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ABSTRACT 

We examine the question of how well the physical properties of clumps in turbulent 
molecular clouds can be determined by measurements of observed clump structures. 
To do this, we compare simulated observations of three-dimensional numerical models 
of isothermal, magnetized, supersonic turbulence to the actual physical structure of 
the models. We begin by determining how changing the parameters of the turbulence 
changes the structure of the simulations. Stronger driving produces greater density fluc- 
tuations, and longer wavelength driving produces larger structures. Magnetic fields have 
a less pronounced effect on structure, and one that is not monotonic with field strength. 
Aligned structures are seen only with low-density tracers, and when the intensity of 
the field is large. Comparing different regions with the same tracers (or conversely, the 
same region with different tracers) can give information about the physical conditions 
of the region. In particular, different density tracers can help determine the size of 
the density fluctuations and thus the strength of the driving. Nevertheless, velocity 
superposition of multiple physical clumps can fully obscure the physical properties of 
those clumps, and short wavelength (compared to the size of the region under analysis) 
driving worsens this effect. We then compare Larson's relationships and mass spectra 
in physical and observational space for the same structure dataset. We confirm previous 
claims that the mean density-size relationship is an observational artifact due to limited 
dynamical range in column density: it is the inevitable consequence presence of a lower 
cutoff in column density. The velocity dispersion-size relationship, on the other hand, is 
reproduced in both physical and observed clumps, although with substantial scatter in 
the derived slope, consistent with observations. Finally, we compute the mass spectra 
for the models and compare them to mass spectra derived from simulated observations 
of the models. We show that, when we look for clumps with high enough resolution, 
they both converge to the same shape. This shape appears to be log-normal, however, 
rather than the power-law function usually used in the literature. 

Subject headings: ISM: clouds, turbulence ISM: kinematics and dynamics, stars: for- 
mation 
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1. Introduction 

This paper tries to understand whether concentrations of emission in observations of position- 
velocity space represent real clumps of gas in molecular clouds; whether the properties inferred for 
the observed objects describe the actual clouds; and what we can say about the mechanisms that 
drive interstellar turbulence by observing molecular clouds with different tracers. 

Several papers have tried to address these points before. Adler &; Roberts (1992) show that, in 
a Galactic disk model, gas complexes in a longitude-velocity diagram do not necessarily correspond 
to real complexes. More recently, based on numerical simulations, Ballesteros-Paredes, Vazquez- 
Semadeni & Scalo (1999) showed that observed clumps frequently come from the superposition 
of several physically disconnected regions in the line of sight at the same radial velocity (see also 
Ostriker et al. (2001),) but not necessarily at the same three-dimensional velocity. Moreover, the 
morphology seen in observed position-position-velocity maps seems to be more representative of 
the spatial distribution of the velocity field in the line of sight than of the distribution of the density 
field, as has been shown in 3D numerical simulations by Pichardo et al (2000). Similarly, Lazarian 
et al (2001) shown that only the velocity field is responsible for the structure existing in the channel 
maps of observational data cubes. 

If observed structures do not match actual structures, it is valid to ask whether statistical 
properties of observed clouds, such as the Larson (1981) relationships, are still valid for actual 
clouds. As Larson (1981) mentions, they may instead be artifacts of the observations. In fact, Scalo 
(1990) argues that the mean density-size relationship is an artifact of the limited dynamical range 
of the observations (see also Kegel (1989).) Vazquez-Semadeni, Ballesteros-Paredes & Rodriguez 
(1997) used two-dimensional simulations of the diffuse interstellar medium to suggest that there 
is no relationship between mean density and size of physical clouds (defined as a connected set of 
pixels in their simulations), but that there is a relationship, although with strong scatter, between 
the velocity dispersion and size, a result confirmed by Ostriker et al. (2001). 

Although as we mentioned, several papers have studied the same questions as the present paper, 
they addressed them only partially, studying for example only objects in physical space (Vazquez- 
Semadeni, Ballesteros-Paredes & Rodriguez 1997), or in observational space (Ostriker et al. 2001); 
or giving physical arguments to explain the observed Larson (1981) relationships (e.g., Larson 
(1981) itself; Kegel (1989); Scalo (1990)). In the present paper we use three-dimensional numerical 
models of molecular clouds to directly compare structures in simulated observations to the actual 
physical structures that produce them. Our goals are: to give observational predictions about the 
nature of the sources of the turbulence; to disentangle what observed properties are representative 
of the true physical entities; and what are more likely to just be the result of superposition of 
multiple structures along the line of sight. Another way of asking this is, what are the differences 
between the observational space and the physical space? 

The plan of the paper is as follows: in §2 we describe the numerical models used here, and 
the methods used to define clumps and to generate simulated observations from model density and 
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velocity cubes. In §3.1 we study the role of the luminosity L, wavenumber k and initial magnetic field 
B in the production of density fluctuations, and give observational criteria to infer the properties of 
the sources of turbulence. Wc discuss the superposition efl"ects, and analyze the relative importance 
of the parameters of the simulations in the superposition effects (§3.2). Section 4 studies Larson's 
relations (mean density vs. size, velocity dispersion vs. size and mass spectrum) and discusses the 
importance of each one of them, by comparing physical and observational space, and §5 draws the 
main conclusions. 



2. Data Analysis 
2.1. Numerical Simulations 



In the present work we analyze driven, supersonic, hydro- and magnetohydrodynamical (MHD) 
models described by Mac Low (1999). Details of the models and forcing can be found there. Here we 
just mention that the models are solutions of the mass, momentum and induction equations using 

ZEUS-3D (Stone & Norman 1992a,b; Stone, Hawley, Gammie, & Balbus 1996) at a resolution of 
128''^ and 256^ zones. They include periodic boundary conditions in each direction and an isothermal 
equation of state, implying that the simulations arc scale free. They are forced at a characteristic 
scale given by the wavenumber k and at a constant energy injection rate L, with an initial magnetic 
field intensity B. 

Table 1 shows the parameters of those simulations. The first column shows the name of the run, 
following the notation in Mac Low (1999)^ In columns 2, 3 and 4 we show the driving luminosity, 
wavenumber and initial intensity of the magnetic field intensities. Column 5 shows the value of the 
root mean square of the density fiuctuations 



ap = 



N 



N-l 



3=0 



1/2 



(1) 



These simulations were performed to model the interior structure of molecular clouds, at scales 
of 0.01-1 pc, densities of n ~ 10^ cm~^, and temperatures of T = 10 K. Under these conditions, the 
isothermal equation of state is a reasonable approximation to the full thermodynamics of the cloud. 
Other works focussing on larger scales (10-1000 pc) have explicitly included heating and cooling, 
a star formation scheme, and so forth (Passot, Vazquez-Semadeni &; Pouquet 1995), but these are 
not necessary to study the scales we are concerned with. A bigger omission in our simulations is 



^The first letter, M or H denotes MHD or hydrodynamical; the second letter gives the strength of the forcing, 
ranging from A through two orders of magnitude in equal logarithmic steps to E; the first number gives the value 
of the driving wavenumber used; and in the magnetic case, the second number shows the ratio va/cs of the Alfven 
velocity to the sound speed, with X representing 10. 
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the lack of self-gravity, which will play a role in the production of dense cores (Klessen, Heitsch & 
Mac Low 2000; Heitsch, Mac Low, & Klessen 2001). Nevertheless, wc find the non-self-gravitating 
simulations useful to study superposition effects, and their consequences in the statistical properties 
observed for clumps. Even more, since our results are consistent with other works that include self- 
gravity, we conclude that its lack in our simulations is not relevant for the conclusions stressed 
here. A similar conclusion has recently been reached by Ossenkopf (2001), who studied wavelet 
transforms of sclf-gravitating MHD turbulence and found that the overall structure was decoupled 
from that of the densest cores. 



2.2. Simulated Observations 

Wc analyze the numerical data first by analyzing the physical variables such as density, veloc- 
ity, and magnetic fields themselves, and second by simulating observations of the models. In pre- 
vious work (Ballesteros-Paredes, Vazquez- Semadeni & Scalo 1999; Ballesteros-Paredes, Hartmann 
&; Vazquez-Semadeni 1999) we represented the line profiles of the simulations as histograms of the 
velocity weighted by density. However, those line profiles did not include the effect of temperature 
in broadening the profiles, and did not account for optical depth effects. Therefore we calculate the 
radiative transfer here, assuming local thermodynamic equilibrium (LTE) for the population of the 
molecular energy levels, which is a sufficiently good approximation to study qualitatively the effects 
of projection and superposition (velocity crowding) of structure in the cloud. In order to model 
differences when observing the same region with different tracers, we took into account that at low 
densities the molecules will be underexcited. To do this we made the very simple assumption that 
the lines arc excited only if the density is above 650 cm~^ for the ^^CO(l-O) transition, 6200 cm~^ 
for 13C0(2-1), and 1.8x10^ for CS(l-O) (Rohlfs k Wilson 1996). In other words, we assume LTE 
if the local density is larger or equal to these density thresholds for each transition. 

Because the simulations are scale-free, the units of density, temperature and length are arbi- 
trary, but the included physics is adequate to represent the behavior of molecular gas with length 
scales between thousands of astronomical units for the densest cores, up to several parsecs. For 
the construction of the ^^CO(l-O), (2-1) and CS(l-O) line profiles, we adopted solar abundances; 
and typical scales for the density and length of n = 1000 cm~^ and L = 0.5 pc, respectively. Since 
the simulations are isothermal, we choose the temperature equal to 10 K. These choices produce 
regions that are mostly optically thin in both ^^CO and CS. 

In the production of line profiles, we fixed the number of velocity channels at 32. This choice 
results in a velocity resolution of 0.06 km s^^ for the less energetic simulations, and of 0.2 km s^^ 
for the more energetic cases. These resolutions are comparable to the velocity resolution of typical 
observations of molecular clouds at FCRAO which is 0.08 km s~^, or in large-scale surveys like the 
full Galactic CO Survey by Dame et al. (2001) which reaches resolutions as high as 0.26 km s"^. 

Since we know the detailed density and velocity distribution, and assume a constant tem- 
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perature, in order to determine the emission from any given molecular transition we just have to 
integrate the corresponding transfer equation, 

^ = n^p{Si, - ly) (2) 

where /j^ is the intensity; is the source function, which we assume to be a blackbody at T = 10 K; 
ds is the distance along the ray; and is the mass absorption coefficient of the line, and has units 
of [cm^ gr^^]. In order to integrate equation (2) along the line of sight we solve, in practice, the 
following equations: 

h = loe-^"'*"' + / e-^^dT^ (3) 

JO 

and 

Kupds (4) 

where /q represents the background intensity and Ti/^tot is the total optical depth of the cloud, 
and p is the local mass density. Since the models have constant elements of volume (zones), the 
integration over s becomes a summation with constant step equal to the zone size. The absorption 
coefficient is 

where A^j is the Einstein coefficient for the i ^ j transition if the density is larger than the threshold 
density for excitation, and zero if the density is smaller; gi and gj are the statistical weights of the 
upper and lower states respectively; rij is the numerical density of molecules in the lower state; and 
^{y) is the line profile, which is assumed to be thermal. These last two parameters take the form 

_-^mol „ / 



-^^{v) (5) 



and 

*^") = ^V^J ^"PV — — J' 

where mo is the mean molecular mass of the gas (assumed to be 2.5 times the mass of the hydrogen 
atom), mmoi is the mass of the emitting molecule, Xmoi is its molecular abundance relative to the 
hydrogen (assumed to be solar), Q is the partition function, and vq is the velocity of the zone, 
which is known. The molecular parameters are given in Table 2 



2.3. Physical and Observational Representation 

In order to represent and compare the physical and observational spaces in two-dimensional 
plots we integrate over one of the dimensions. In particular, when wc integrate in the x direction 
the physical space, the corresponding observational cube is integrated in the z direction. This 



-6- 



will allow us to compare, whether a clump in a position-velocity diagram has a counterpart in the 
position-position diagram. Then, when we label a plot as z — x{l — 0) or z — x{2 — 1), we show the 
column density field integrated in the y direction, after zeroing all those points where the density 
is smaller than the density threshold (see Table 2.) The corresponding observational plot will be 
z — Vx{l — 0) OT z — Vx{2 — I) respectively, and then represents the simulated observation of the 
z — X map as seen by an observer located at the left of the plot, assuming ^^CO(l-O) or ^^C0(2-l) 
emission as described above. 

2.4. Defining a Cloud 

The first problem that appears when we study the properties of clouds and clumps is how to 
define them. For instance, if we are interested in studying large complexes and their substructure, 
we define a complex as a connected set of points with intensity above an arbitrary threshold (Dame 
et al 1986; Vazquez- Semadeni, Ballesteros-Paredes &; Rodriguez 1997). When we work in physical 
space, the "intensity" is the actual density. If we work in the observational space, the "intensity" is 
the brightness temperature. If instead, we are interested in studying individual clumps, we define a 
clump as a connected set of points below a local maximum following the intensity only downwards 
until the threshold is reached. This scheme has been implemented by Williams et al. (1994) in 
a code called clumpfind. A single large cloud with internal structure in the first scheme may 
be counted as several smaller clumps in the second scheme. Other schemes are also used. For 
example, one can assume that clumps have a Gaussian distribution of intensity (Stutzki & Giinsten 
1990), so that a clump defined as single by clvimpf ind can be decomposed into several Gaussian 
substructures. 

In the present work we will use clumpf indfor simplicity. Differences between methods will be 
discussed in future contributions, but it is worth noting that, as Sheth et al (2000) has shown, 
these methods can not identify clumps larger than a few resolution elements. 

3. Results 

3.1. Role of the parameters 

We analyze how the parameters of the simulations (driving strength and wavelength, magnetic 
field strength) affect the structure in both observational and physical space. In particular, we 
attempt to give observational criteria to help determine the equivalent parameters of the physical 
mechanism driving the observed turbulence. 

We begin by comparing models with the same driving wavcnumber and magnetic field, but 
different driving strengths. In Table 1 we see that the larger the driving strength L, the larger 
the density fiuctuations ap. This occurs because the rms Mach number in the simulations depends 
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on the driving strength and the wavenumber as M oc (L/k)^/^ (Mac Low 1999), and the density 
contrast across an isothermal shock depends on the Mach number as Sp/p oc M^. 

Figure 1 shows maps of both the physical column density and the observed emission for two 
hydrodynamical models with the same wavenumber (k = 8), but driving strength different by a 
factor of 100: model HAS (L = 0.1) and model HE8 (L = 10). The main difference between these 
extreme cases is that HAS has larger column densities of low-density tracers, and lower column 
densities of high-density tracers. Another clear difference is that in high-density tracers HAS shows 
well-separated, roundish cores, while HES shows a more filamentary structure with few compact 
cores. 

These differences can be understood as follows. When the turbulence is weakly driven (HAS), 
density fluctuations are small. There will be few regions with density below the density threshold 
for the low density tracer (^^CO(l-O)), so most of the cube will emit in that line, but conversely 
only a few places will have densities higher than the threshold for the high-density tracer (^^C0(2- 
1)), so only a small column density is seen in that line. Conversely, for strong driving (HES), 
the fluctuations around the mean are larger, and more regions will have densities either below the 
threshold for ^^00(1-0) or above the threshold for ^^00(2-1), so that the mean column density for 
^^CO(l-O) is lower, but the mean column density for ^^C0(2-l) is larger. Stronger driving leads 
not just to higher densities, but more extreme fluctuations, as indicated by a greater fraction of 
material at both high and low densities. 

Since structures in the simulations are produced by the convergence of turbulent flows, the 
smaller the driving wavenumber, the larger the dominant structures. We can sec this effect in 
Figure 2, where we show four maps for HC2 and four maps for HCS, both with the same driving 
strength (L = 1), but different wavenumber (A; = 2 and S respectively). From these plots we can 
see that the density structure is more homogeneously distributed in the case of high wavenumbers 
(HCS), in both physical and observational spaces. On larger scales, this suggests that large molecu- 
lar clouds might be the result of turbulence driven at large scales in the diffuse interstellar gas (see, 
e.g., Ballesteros-Paredes, Hartmann & Vazquez- Semadeni (1999); Hartmann, Ballesteros-Paredes 
& Bergin (2001)) , and not the agglomeration of small clouds by the action of the self-gravity, as 
in some models of molecular clouds (e.g., Scoville & Hersh (1979); Kwan (1979); Norman & Silk 
(19S0)) 

The effect of magnetic fields in our models is consistent with previous work using simulations of 
the diffuse interstellar medium by Passot, Vazquez-Semadeni &; Pouquet (1995), who showed that 
weak magnetic fields decrease the size of the density fluctuations compared to the hydrodynamic 
case. For intermediate magnetic fields, the density fluctuations may increase above the hydrody- 
namic and weak field cases, and finally, for very strong magnetic fields, the density fluctuations 
decrease again. (The actual diagnostic they used, as shown in their Fig. S, is the star formation rate 
as a function of the magnetic field strength; but in their scheme the star formation rate depends 
directly on the size of the density fluctuations). 
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In Figure 3 we demonstrate this behavior. The simulations used have intermediate driving 
strength (L = 1) and wavenumber (k = 4), and three different initial magnetic field strengths of 
B = (HC4), B = 0.1 (MC41), B = 0.5 (MC45). Wc can understand the dependence of the density 
fluctuations on the magnetic field strength as follows: weak magnetic fields are strongly tangled 
by the flow, so that they have a nearly isotropic magnetic pressure that prevents compressions and 
large density fluctuations, while stronger fields have an anisotropic magnetic pressure that allows 
larger compressions along the mean field lines, and thus larger density fluctuations. 

Magnetic fields do not noticeably affect the structure of the maps unless they are strong 
enough to introduce some anisotropy. For example, in Figure 4 we show -"^^ CO (1-0) maps of three 
different simulations, integrated parallel (z), and perpendicular (x) to the initial field direction. 
From bottom to top, the relative importance of the magnetic field compared to the luminosity of 
the simulation increases. We show: ME21 (large driving strength, small wavenumber and small 
magnetic field); MA81 (low driving strength, large wavenumber and small magnetic field); and 
MA4X (small driving strength, intermediate wavenumber and large magnetic field). Only the 
strong field case MA4X (top panels) shows a morphology different than the hydrodynamic case. 
In this case, a marked anisotropy is visible, with the projection perpendicular to the magnetic 
field (upper left panel), showing low-density structures aligned mainly in the direction of the field 
lines while the projection along the field line (upper right panel) presents regular structures, with 
no preferred direction. Interestingly, in a higher density tracer, there is no clear morphological 
difference between projections along or perpendicular to the initial field (see Fig 5.) 

The maximum values of the column density are similar in both physical and observational 
projections for small magnetic field intensity. Nevertheless, in the strong-field case (MA4X, see 
upper panels in Fig. 4 and 5) the column density integrated perpendicular to the mean field (right 
panels) is higher than the column density integrated parallel to it (left panels), as can be seen 
in the grayscale bars. This effect can be explained in the same way as that shown in Figure 3: 
strong magnetic fields only allow compressions parallel to the field. Therefore, when we observe 
along the magnetic field lines we do not see high column density contrasts, while when we observe 
perpendicular to the field lines, we see higher contrasts in the compressions along the field lines. 
As a summary, strong fields can produce anisotropy and compressed structures perpendicular to the 
field lines, but weaker fields merely reduce the maximum density contrast somewhat. 

3.2. Superposition Effects 

Velocity crowding contributes substantially to the generation of clumps in observational space, 
so observed clumps frequently contain emission from physically separated regions (Ballesteros- 
Paredes, Vazquez-Semadeni & Scalo 1999; Pichardo et al 2000; Ostriker et al. 2001; Lazarian et 
al 2001). We demonstrate this effect using a typical MHD simulation: MC41. We choose this run 
because it has intermediate values of driving strength {L = 1) and wavenumber {k = A), and it 
clearly shows the effects of superposition. 
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Figure 6 shows physical and observed maps for ^^CO(l-O) and (2-1). We see that clumps 
in real space (letters A, B, and C in panel y — .x(l — 0)) do not necessarily have a counterpart 
in observational space. Clumps in observational space, on the other hand, (letter D in panel 
— a;(l — 0)) do not necessarily come from isolated regions in real space, but have contributions 
from many different regions along the same line of sight. In Figure 6 we plot dotted lines that 
show the places where the emission of a physical clump lies in observational space, and where the 
emission from an observational clump is generated in physical space. For reference, we use the 
same lines in the ^^CO(l-O) map as in the ^^C0(2-l) map. 

For the same data set we made another set of maps using higher density tracers. Figure 7 
shows the emission of CS(l-O). The same lines are shown as in Figure 6, to indicate the location of 
the large CO (1-0) core. In this case, the illusory observational clump D has almost disappeared, 
and the physical clump A is observed more clearly in its observational counterpart. We also note 
that the position that we would attribute to the cores in different tracers is slightly different. 

In order to quantify the effects of superposition, we now compute the typical number A^ciump 
of physically disconnected regions that contribute to the emission of a single observed clump. To 
do so, we developed an algorithm that brackets each observed clump in position and velocity and 
calculates the number of connected regions in physical space within the bracketed region in phase 
space. 

In Figure 8 a we show the influence of driving strength and wavenumber on the amount of 
superposition. We find that high-density tracers behave differently than low-density tracers. In 
the case of a high-density tracer such as CS, A^^^iump 2-3, and is almost constant for different 
driving strengths and wavenumbers. On the other hand, in the case of a low-density tracer such as 
^^CO, which emits from most of the cube, superposition increases with driving wavenumber. This 
is a natural consequence of the way we forced the models, since short-wavelength driving produces 
multiple physical regions at the same velocity. There is also a marginal tendency for strongly-driven 
models to exhibit more superposition than weakly-driven ones. This occurs because the greater 
velocity dispersion in the strongly-driven model produces clumps with greater internal velocity 
dispersion that more easily overlap. We conclude that higher-density tracers will better disentangle 
velocity crowding, as is shown by the results using CS(l-O), though still by no means completely. 

In Figure 86 we show the relative influence of the magnetic field by examining A'ciump as a 
function of the ratio of the driving strength to the initial magnetic field strength L/B. We find that 
the amount of superposition is a bit larger when we observe along the mean field direction (z) than 
when we observe perpendicular to it. The effect is small, but present, and it is a consequence of 
the greater compressions parallel to the field lines producing more condensed clumps with smaller 
cross-sections for superposition when observed perpendicular to the field lines. 

A final word concerning superposition is in order: although it is true that high spatial resolution 
helps to segregate large clouds into their clumps, high velocity resolution does not necessarily help 
to disentangle superposition. It only helps if there is a monotonic velocity gradient along the line 
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of sight, such that two clumps with similar velocities merge at low velocity resolution, and can be 

separated at higher velocity resolution. If clumps are produced in regions of converging flow in 
a larger turbulent region, however, as is likely the case in real clouds, then the velocity gradients 
will not be monotonic, and physically disconnected regions will often share the same line-of-sight 
velocity, so increasing velocity resolution will not help. 

As a summary, we conclude that velocity crowding is likely to always occur, implying that 
cores in observational space are not single entities in real space, and that it is larger for low-density 
tracers, large driving wavenumbers, and when we look along the mean field lines. 



4. Larson's Relations in Numerical Simulations 

If observed clumps have contributions from multiple physical regions, it is important to under- 
stand whether the relationships reported for observed clumps in molecular clouds have anything to 
do with the properties of the actual physical clumps. 

Larson (1981) studied the dependence on size of the mean density, velocity dispersion, and 
mass spectrum of the clouds in a sample of observational data taken from the literature. He found 



{p) oc i?", (8) 
5v oc i?"*, (9) 



and 

dN 



d\ogM 



oc M^. (10) 



The most commonly quoted values in the literature are a ^ —1, m ~ 1/2, and 7 ~ — 0.5. In 
particular, values of ck = — 1 and m = 1/2 have been attributed to virial equilibrium^ (Larson 1981; 
Myers &; Goodman 1988; Vazquez-Semadeni, et al 2001). However, there is some discrepancy in 
the values reported (Carr 1987; Loren 1989), especially for high latitude clouds, where the velocity 
dispersion-size relationship has been questioned (LaRosa, Magnani k. Shore 1999). 

Larson (1981) himself already suggested that the relationships he identified might be due to the 
limitations of the observations, but it was Kegel (1989) who first demonstrated that the observed 



^What is often called virial equilibrium should actually be called energy equipartition, since (a) observations are 
not able to measure all the terms in the virial equation, such as the surface and time-derivative terms, to prove 
virial equilibrium, and (b) it has been shown that clouds in (two-dimensional) numerical simulations of the diffuse 
interstellar medium are in energy equipartition (Ballesteros-Paredes & Vazquez-Semadeni 1995) but not necessarily 
in virial equilibrium (Ballesteros-Paredes & Vazquez-Semadeni 1997). 
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mean density-size relationship could be due to observational effects, and that the observed and 
physical properties such as radius or volumetric density may be quite different. Scalo (1990) 
showed that CO and extinction saturate at roughly the same column density, which forms the 
upper envelope of the mean density-size relationship. He also discussed how the surveys are biased 
due to the limitations of instrumental sensitivity and saturation. Vazquez- Semadeni, Ballesteros- 
Paredes & Rodriguez (1997) reproduced the mass spectrum, as well as the velocity dispersion- 
size relationship (although with large scatter) for clumps in physical space, using two-dimensional 
models. They reported the lack of a mean density-size relationship, confirming numerically the 
analysis by Kegel (1989, sec also the discussions in Larson (1981) and Scalo (1990)) in the sense that 
there are clouds with small sizes and low column density that will be undetected in observational 
surveys. More recently, Ostriker et al. (2001) reported a slope close to —1 for the mean density- 
size relationship derived from simulated observations of their models^, and a very flat correlation 
between velocity dispersion and size (m ~ 0.1). As they mentioned, both results are a consequence 
of their clump finding method (based on what they call regions of contrast), which samples the 
entire line of sight. 

None of these studies attempting to explain the origin of the Larson (1981) relationships have 
compared simulated observations to the actual physical density and velocity distributions for the 
same numerical model. We now do that, using the three-dimensional simulations of molecular 
clouds that we have described above. 

We must first define the size of a clump, both in observational and in physical space. For 
clumps drawn from simulated observations we take the circular radius, defined as Rdrc = (^/^r)^^^, 
where A is the projected area of the clump on the plane of the sky, for comparison with WGB94. For 
physical clumps we take the geometric radius, defined as R= (i?^ + Ry + i?^)^/^. Other definitions 

are equally valid, but can lead to differences of a factor up to half an order of magnitude when 
comparing different size definitions, and differences of as much as a factor of 100 when calculating 
mean densities by dividing the total mass by the size^. 

Finally, since the results discussed in the present section do not appear to depend on the 
parameters of the runs we analyze or (in the magnetic case) on the projection with respect to the 
mean magnetic field, we only present results for run HC8-256 (see Table 1), an available high- 
resolution run (256^ zones). We search for clumps in physical space using clumpf ind, considering 
only regions with a density larger than 650 cm~^, the critical density for CO (1-0) emission. The 
results do seem to depend on the clumpfinding scheme used and on its parameters. For example, 
the number and size of the clumps found in clumpf ind depends on the level of refinement used in 
clumpf ind: at high refinement, large clumps can be split into multiple smaller clumps, changing not 
only the slope of the mass spectrum, but also the dynamical range of size in the scaling relations. 



^Actually, they report a slope of 2 in the mass-size relationship, but the interpretation and result are the same. 

^Assuming LTE. Non-LTE effects will also increase the errors (Padoan, Juvela, Bally, & Nordlund 2000; Ossenkopf 
2001). 
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4.1. Mean Density-Size relationship 

In Fig. 9 (upper panel), we show the mean density-size relationship for clumps in three- 
dimensional physical space. We note three points about this Figure. First, the relation between 
mean density and size given by equation (8) is not seen at all, confirming the results of VBR97 
from two-dimensional models. Second, there is a minimum density below which there are no clumps 
identified. This minimum is just given by the density threshold we used in clumpf ind. Third, even 
though the simulations exhibit a large dynamical range in density (pmax/Pmin ^ 3.5 x 10^), the 
dynamical range in the mean density-size relationship is small, because in constructing such a plot 
we choose the clumps around the local density maxima. 

In the middle and lower panels of Fig. 9 we show the mean density-size relationship for clumps 
in simulated observational maps of the same run, integrated along the x axis. The middle panel 
is the one obtained by using CS(l-O), and the lower panel is the relationship obtained by using 
^^CO(l-O). The observed clumps do exhibit approximately the relationship given by Larson (1981), 
despite the lack of correlation exhibited by the physical clumps in this model. We conclude from 
this demonstration that the observed density-size relationship (eq. [8]) is an observational artifact. 

Two mechanisms have been suggested to explain density-size relationship with slope ~ — 1 in 
observations. First, it might be due to the selection effect described by Kegel (1989): only clouds 
with intensity exceeding the noise threshold determined by the instrumentation will be detected, 
effectively setting a column density cutoff, rather than the physical density cutoff imposed in the 
physical density-size relationship. A constant column density cutoff produces a cutoff with slope 
— 1 in the mean density-size plane (middle and lower panels), just as the constant physical density 
cutoff in physical space produces a flat cutoff. 

Second, it might be due to the limited dynamical range of the observations (Scalo 1990). To 
understand this, lets consider the simulations: if we had clumps with constant column density N, 
we would infer a mean density {p{R)) = NR^^, giving exactly a slope of —1. In the simulation 
analyzed here we do not have a constant column density, but the column density varies by only a 
factor of 30 (see Burkert &; Mac Low 2001), far less than the variation in physical density. Even 
worse, from this narrow distribution, the selected clumps have an even narrower column density 
distribution: on the one hand, as we already mentioned, there is a column density cutoff given by 
the minimum intensity used in clumpf ind; on the other hand, even the emission of the brightest 
clump does not come from the whole line of sight, and its column density is not the largest column 
density available in the simulation. The end result is that selected clumps naturally end up with 
column density constant to within an order of magnitude, consistent with the observed scatter in 
the middle and lower panels of Fig. 9. 

Similar arguments can be applied to the observational data. For example, when observing 
emission lines from molecules, there is a limited dynamical range in density available in which 
any particular tracer can be excited: it will not be excited at too low density, and it will become 
optically thick at too high a density Scalo (1990). At IRAS wavelengths. Wood, Myers & Daugherty 
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(1994) claim constant column densities in IRAS observations of dust, where the dynamical range of 
the data should have allowed detection of any significant variations. However, Vazquez- Semadeni, 
Ballesteros-Paredes &; Rodriguez (1997) showed that for IRAS observations, only the external layers 
of the clouds are heated by the diffuse UV radiation, so that the emission at 60 and 100 /im actually 
comes from surface layers that are indeed expected to have nearly constant column density. 

Finally, we point out that deviations from the mean density-size relationship have also been 
found in observations (Carr 1987; Loren 1989; Falgarone, Puget & Perault 1992), suggesting that 
with sufficient care, the trap of apparently constant column density can be overcome. 

4.2. Velocity Dispersion-Size Relationship 

In Figure 10 we show the velocity dispersion-size relationship for clumps in physical and obser- 
vational space. In the case of physical space (upper panel) the typical velocity dispersion is of order 
of 0.2 km s~^, a value that coincides with the value of the fluctuations of the overall velocity field 
for the run analyzed here (HC8-256). The slope fit to the data is 0.3 it 0.38, smaller than the value 
of 1/2 suggested for compressible turbulent behavior (see, e.g., VBR97), but with a strong scatter. 
In fact, in order to reduce substantially the uncertainty in the determination of the slope, we should 
have at least two orders of magnitude of dynamical range in sizes if the scatter is about one order 
of magnitude, as it is the case. But with clumpf ind, it is difficult to have a large dynamical range 
(see, e.g., upper right panel of Fig 7 in Williams et al. (1994)), since a reasonable refinement of the 
intensity contouring makes large clumps split up into several pieces (Sheth et al 2000). 

On the other hand, when we analyze observational data as in Fig 10 (lower panel), additional 
problems play a role. We now have the possibility that substantial superposition is occurring, in 
particular for low-density tracers. If this is so, the velocity dispersion of a single observed clump 
will reflect the velocity dispersion of several, physically disconnected regions in the same line of 
sight. This will also help to produce both large mean values of Sv (0.8 km s~^, compared to the 
0.2 km corresponding to the mean velocity fluctuations in the simulation) and low values of 
the slope in the velocity dispersion-size relationship (Fig. 10, lower panel). 

In contrast, for high density tracers (CS(l-O), Fig.lO, middle panel), we obtain a steeper 
relationship (0.46), and the mean velocity dispersion drops back to ~ 0.2 km s~^. Three reasons 
why this occurs are: in the case of high density tracers there is weaker superposition as shown 
in Fig. 8; clumps are substantially smaller than in the low-density case, giving smaller velocity 
dispersion; and the sizes span about one and a half orders of magnitude, favoring disentanglement 
between the scatter and the actual trend. 
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4.3. Mass Spectrum 

Now we turn to the clump mass spectrum. Before showing the mass spectra for the simulations, 
we have to think about the nature of the turbulent density structure and how clumpf ind (or any 
other clumpfinding algorithm) works. For instance, in a driven turbulent medium, such as the 
simulations we examine here, we expect density fluctuations at all scales above the dissipation 

scale. It is not at all clear when to stop counting substructure inside structure. In other words, 
should we look for all the stones in the small hills in every single mountain of the mountain chain, 
or only count the larger mountains? If only the larger mountains, how do we define a lower limit? 

Intuitively, we might think that the dependence of clump numbers on clumpf ind parameters 
would be reduced in regions where clumps are well defined and reasonably isolated, as in an isolated 
star-formation region where cores have suffered some amount of collapse. To show that this is not 
the case, we show in Figure 11 a contour map of a slice through the center of model HC8. We 
use five and ten isocontours between an arbitrary density threshold of Pmax/S and the maximum 
density Pmax- As we can see, the density structures seem to be well-separated and reasonably 
well defined. In Figure 12 we present an enlargement of the lower-right corner of Figure 11, again 
using five and ten isocontours. We see that, even for the same density threshold, by changing the 
number of isocontours we find five more peaks (see right panel in Fig. 12) that will result in five 
new clumps in clumpf ind. This happens in just a small region containing what we would think 
of as well-defined structures, at a single slice in z. But in practice, we have 3D structures, so that 
every structure is far more complicated than apparent at first glance. Even without considering 
whether observational and physical clumps are the same or not, the mass spectrum appears to be 
more an artifact of the manner in which we count clumps than an actual physical characterization 
of the structure, unless we count them exhaustively all the way down to the numerical or physical 
dissipation scales below which there is no longer well-defined structure. We note also that in the 
case of self-gravitating simulations the result is the same (Osscnkopf et al. 2001), and that similar 
results will hold for observations: at low resolution the structures look soft, but once we increase 
the resolution and the sensitivity, structure emerges inside the structure, as is clear from Figure 4 
in Dame et al. (2001). 

Preliminary work by Nordlund (2002) shows that, by setting the density threshold in physical 
space to the average density, a large enough number of contour levels gives a log-normal mass 
spectrum, independently of the number of contours used. We compute the mass distribution using 
this prescription for run HC8 smoothed at a resolution of 64^ for physical space, and 64^ x 32 
velocity channels for observational space in order to avoid excessive computation time and memory 
usage by clumpf ind. The results are shown in Fig. 13. The upper panel corresponds to physical 
space, and the lower panel corresponds to observational space, where instead of density, we used 
intensity contouring. The results are clear: for large number of isocontours the shape of the mass 
spectrum converges. Good agreement is observed in the large-mass end of the distribution, although 
slower convergence is observed at the low-mass end. Convergence to a factor of two for the full 
spectrum appears to be reached for 64 or more contours. 
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In Fig. 14 we compare the mass spectrum calculated in physical space (bold line) to that from 
observational space (thin line) using 512 isocontours. The dotted parabola shows a log- normal fit 
that reproduces fairly well the large-mass end of the spectra. Such log-normal behavior of the 
observed spectrum might explain observations obtained for dense cores in the last few years. In 
particular, Testi & Sargent (1998) found a power-law mass spectrum with slope of —2.1 for star- 
forming cores in Serpens, substantially larger than the standard —1.5 value. Moreover, Motte, 
Andre &; Neri (1998) found a change in slope at low mass: the less massive cores exhibit slope of 
— 1.5, considerably shallower than the —2.5 slope found for the more massive cores. Finally, we note 
that in these works, the mass spectrum was obtained using dust continuum observations with high 
signal-to-noise ratios. For present-day molecular-line observations, typical signal-to- noise ratios of 
10 or less do not allow small contour separation, since contouring at less than the signal-to-noise 
ratio is meaningless (Williams et al. 1994). Mass spectra derived from such observations are likely 
to be far from physical (Fig. 13). 

5. Summary 

We analyze a set of three-dimensional MHD numerical simulations at intermediate (128^ zones) 
and higher (256"^ zones) resolution described by Mac Low (1999), in which the parameters driving 
strength (L), driving wavenumber (k), and initial magnetic field (B) are varied. We show that: 

• Density fluctuations are primarily controlled by the driving luminosity: the larger the lumi- 
nosity, the larger the fluctuations seen in the maps. Morphology of the density structure, 
on the other hand, is primarily controlled by the driving wavelength: the larger the driving 
wavelength, the larger the structures seen in the maps, and vice- versa. 

• The magnetic field only produces detectable effects in our simulated molecular line obser- 
vations if the fields are strong and in the plane of the sky, and the structure is observed in 
low density tracers. When observing magnetized clouds, it is slightly easier to obtain larger 
column densities if the field is in the plane of the sky than if it is perpendicular to it. 

• Simulated observational maps of the models show strong superposition of structures along 
the line of sight. This implies that single clumps in observational space arc frequently the 
superposition of multiple clumps in physical space. The apparent properties of observed 
clumps are poorly related to the physical properties of their constituent objects if strong 
superposition occurs. 

• We explore the effect of the simulation parameters on the generation of strong velocity crowd- 
ing, finding that large driving strength and large wavenumber tend to favor velocity crowding. 
So far, the only more-or-less reliable method to ensure there is no substantial superposition 
of clumps is to use high-density (close to the maximum density) tracers. Nevertheless, even 
in those situations superposition can occur. 
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• Larson's (1981) mean-density to size relationship does not exist in physical space, but does 
occur in our simulated observations. It appears to be an observational artifact, whose expla- 
nation we discuss. 

• Larson's velocity-dispersion to size relationship is reproduced in both physical clumps and 
simulated observations, though with substantial scatter in the slope. More dynamical range in 
the simulations would probably help, although our clumps already span an order of magnitude 
in size with half an order of magnitude in scatter. 

• The mass spectra converge to a particular shape only when we use large number of contours 
(above 32 for the high mass end, but above 64 for the whole spectrum), corresponding to 
an extraordinarily large signal-to-noise ratio in the observations. Rather than a single power 
law, they seem to follow a log-normal distribution. 
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Fig. 1. — Physical {z—x) and observational (z—Vx) maps for two models with the same wavenumber 
{k = 8) and no magnetic field {B = 0), but different values of the driving strength: (a) HAS 
(L = 0.1), and (b) HE8 {L = 10). For weaker driving, the density fiuctuations axe small and the 

entire cube contributes to the emission in low density tracers (HAS z — a; (1-0) and z — Vx{l-0)), but 
only a few places contribute to the emission in high density tracers (HAS z — x{2-l) and z — Vx{2-l)). 
In contrast, for stronger driving, pixels with density below the threshold do not contribute, giving 
comparatively lower column density values in low density tracers (HE8 z — x(l-O) and z — 14(1-0)), 
but substantially more regions are at high densities, giving higher column densities in high density 
tracers (HE8 z - x{2-l) and z - T4(2-l)). 

Fig. 2. — Physical and observational maps for two runs with the same driving strength (L = 1) 

and no magnetic fields (S = 0), but large and small values of the wavenumber: (a) HC2 {k = 2), 
and (b) HC8 {k = 8). The larger the wavenumber, the smaller the size of the dominant scale for 
density structures in the simulation. 

Fig. 3. — Root mean square of the density fiuctuations (dp) as a function of the initial magnetic 
field strength of the simulations. Low magnetic field strength prevents density fiuctuations more 
efficiently than either zero or larger-strength fields. 

Fig. 4. — Physical {z — x(l-O) and y — x(l-O)) maps for large-scale, large intensity driving and 
small magnetic fields (ME21, lower panels); small-scale, small intensity driving and small magnetic 
fields (MA81, middle panels); and intermediate scale, small intensity driving and strong magnetic 
fields (MA4X, upper panels) integrated along the y (left panels) or z (right panels), where the 
initial magnetic field runs parallel to the z-axis. The relative importance of the magnetic field 
(compared to the luminosity) increases upwards. Column densities are similar in both projections 
(perpendicular and parallel to the mean magnetic field lines), except if the magnetic field is strong 
(upper panels). In this case, projections along the magnetic field gives smaller column density 
than the projection perpendicular to it. As explained in the text, strong magnetic fields allow only 
compressions parallel to the field, and then only when observing perpendicular to the field lines we 
get high density contrasts. 

Fig. 5. — Same as Figure 4 but for (2-1). Note that the total column density is still typically larger 
when we integrate perpendicular to the field lines than along them. Nevertheless, differences in the 
morphology are more subtle, showing that the alignment of the structures to the magnetic field 
occurs preferentially at low densities. 

Fig. 6. — Physical and velocity space maps for ^^CO(l-O) and (2-1) in run MC41 (L = 1, A; = 4 
and B = 0.1). Clumps in physical space (A, B, C) do not necessarily correspond to clumps in 
the observed velocity space. Observed clumps in velocity space (D) are not necessarily formed by 
emission from a single region in physical space. 
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Fig. 7. — Physical and velocity space maps for the CS emission of run MC41. Since CS traces 
higher densities, the clumps are more separated than in the ^^CO emission shown in Figure 6. 
Note that the position of the cores A, B, C, D has changed slightly respect to their position in 
Fig. 6 

f8b.eps 

Fig. 8. — (a) Mean number of physical regions contributing to each clump in ^^CO(l-O) and in 
CS(l-O) (A^ciump) versus log2 k. The different types of points denote different driving luminosities 
L. Superposition accounts better for large wavenumbers, low-density tracers, strongly driven sim- 
ulations, (b) A'^ciump for models both parallel (z) and perpendicular (x) to the mean magnetic field 
direction vs ratio of driving strength to initial field strength L/B for models MA4X (weak driving, 
strong field), MC4X (stronger driving, strong field), MC41 (stronger driving, weak field). Error 
bars show the standard deviation. 

Fig. 9. — Mean density-size relationship for physical clumps in physical (upper panel); and simu- 
lated observational clumps in observational coordinates (middle and lower panels) . The dotted line 
has a slope of a = — 1. In physical space we find no correlation, verifying the results by VBR97, 
but nevertheless the simulated observations show such a correlation, as found by Larson (1981) and 
many others. The selection of two different density tracers was chosen to show that the apparent 
correlation does not depend on the selection of the density threshold. 

Fig. 10. — Velocity dispersion-size relationships for run HC8-256 in physical space (upper panel), 
and observational space (middle and lower panels) using CS(l-O) and ^^CO(l-O) respectively. The 
dotted line has a slope of 1/2, the expected value for a turbulent medium dominated by shocks 
(see, e.g., Vazquez-Semadeni et al. 2000). The solid line is the least squares fit to the data points, 
with a slope m and its uncertainty shown in each frame. 

Fig. 11. — Isocontours for the density field p > Pmax/5 at z = 64 pixels using (a) five isocontours, 
and (b) ten isocontours. Note that the clumps seem to be well-defined and isolated. 

Fig. 12. — Enlargement of the lower-right corner of Fig. 11. We use again (a) five, and (b) ten 
isocontours. Note that many new clumps appear with the increased number of contours, even for 
well defined, isolated clumps. 

Fig. 13. — Mass distribution of physical (upper panel) and observational (lower panel) clumps in 
run HC8, computed using clumpf ind with density threshold equal to the mean density (physical 
space) or mean intensity (observational space), and varying the number of contours: 8 (dotted line), 
32 (thin, solid line), 128 (medium bold solid line), and 512 (bold solid line). After ~ 64 contours 
(not shown here for clarity purposes) , the shape of the whole mass spectrum converges reasonably 
well. Below that resolution, the final number of clumps, and especially their distribution for the 
low-mass end varies strongly with the contouring chosen. 
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Fig. 14. — Mass distribution (solid lines) for physical (thin) and observational (bold) spaces using 
clumpf ind with 512 contours and a threshold equals to the mean density or intensity, respectively. 
Note that the mass spectrum inferred for observations and simulations converges to the same shape, 
a log-normal distribution, as shown by the dashed-line parabola. 



Table 1. Properties of the simulations. 



Name 


T 

L/ 


K 


r> 
JD 




MA81 


0.1 


8 


0.1 


0.68 


MC81 


1 


8 


0.1 


1.21 


MC85 


1 


8 


0.5 


1.56 


MC45 


1 


4 


0.5 


1.83 


MA4X 


0.1 


4 


1 


1.27 


MC4X 


1 


4 


1 


2.21 


MC41 


1 


4 


0.1 


1.35 


ME21 


10 


2 


0.1 


2.32 


HC4 


1 


4 





1.97 


HC8 


1 


8 





1.40 


HC8-256 


1 


8 





2.00 


HAS 


0.1 


8 





.71 


HB8 


0.3 


8 





1.05 


HD8 


3 


8 





1.82 


HE8 


10 


8 





1.97 


HE4 


10 


4 





2.39 


HC2 


1 


2 





2.21 


HE2 


10 


2 





2.50 



Note. — All runs are at a res- 
olution of 128^ except run HC8- 
256, which is the equivalent to 
run HC8 but at a resolution of 
256^ pixels 
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Table 2. 



Molecule 


Frequency 




Dipole Moment 


^mol 


Threshold density 


References 




(s-i) 


(s-i) 


(Debye) 


(relative to H2) 


for excitation 




13CO{1-0) 


1.1 X 10" 


7.6 X 10-8 


0.112 


1.43 xlO-« 


650 cm-3 


1,2 


i3CO(2-l) 


2.2 X 10" 


6.5 X 10-'^ 


0.112 


1.43x10-6 


6200 cm-3 


1,2 


CS(l-O) 


4.9 X 10^° 


2.1 X 10-8 


1.957 


1.0 xlO-i° 


18000 cm-3 


1,2 



iRohlfs & Wilson (1996) 
2 Gomez & D'Alessio (2000) 
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